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fidelity and double orthogonality functions are compared as methods for the determination of the 
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using complex scaling formalism. 
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I. INTRODUCTION 



In the past few years the apphcation of quantum information concepts to some long 
standing problems has led to a deeper understanding of those problems jl| and, as a conse- 
quence, to the formulation of new methods to solve (or calculate) them. For example, the 
simulability of many body problems is determined by the amount of entanglement shared 
between the spins of the system 

There is a number of quantities that can be calculated in order to analyze the information 
carried by a given state including the entanglement of formation [sl, the fidelity {4] and 
several kinds of entropies, entanglement witnesses and so on. Of course it depends on the 
problem which quantity is more adequate, or accessible, to be calculated. 

In the case of atomic or few body systems with continuous degrees of freedom a rather 
natural quantity is the von Neumann entropy, it has been used to study a number of prob- 



lems: Helium like atom 



a. 



7|, the dynamical 



6| , generation of entanglement via scattering 
entanglement of small molecules js], and entanglement in Hooke's atom [ol 

In quantum dots, most quantum information studies focus in the amount of entanglement 



carried by its eigenstates 
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ll| , or in the controllability of the system [12] . Both approaches 



are driven by the possible use of a quantum dot as the physical realization of a qubit[13|. 
The controllability of the system is usually investigated (or performed) between the states 
with the lowest lying eigenenergies fl^ . 

Besides the possible use of quantum dots as quantum information devices there are pro- 
posals to use them as photodetectors. The proposal is based on the use of resonance states 
because of its properties, in particular its large scattering section compared to the scattering 



sections of bounded states 



16|. 



The resonance states are slowly decaying scattering states characterized by a large but 
finite lifetime. Resonances are also signaled by sharp, Lorentzian-type peaks in the scattering 
matrix. In many cases of interest where complex scaling (analytic dilatation) techniques 
can be applied, resonances energies show up as isolated complex eigenvalues of the rotated 
Hamiltonian [l5|. Under this transformation the bound states remain exactly preserved and 
the resonance states are exposed as functions of the ^rotated Hamiltonian. Resonance 
states can be observed in two-electron quantum dots 

Recently, a work by Ferron, Osenda, and Serra 



16 



171 ] and two-electron atoms 



19| studied the behavior of the von 
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Neumann entropy associated with C? approximations of resonance states of two electron 
quantum dots. In particular that work was focused in the resonance state that arises when 
the ground state loses its stability, i. e. the quantum dot does not have two electron bounded 
states any more. Varying the parameters of the quantum dot allows the energy to cross the 
ionization threshold that separates the region where the two electron ground state is stable 
from the region where the quantum dot loses one electron. 



In reference 



19| it was found that the von Neumann entropy provide a way to obtain the 



real part of the energy of the resonance, in other words, the von Neumann entropy provides 
a stabilization method. The numerical approximation used in 19| allowed to obtain only a 
reduced number of energy levels (in a region where the spectrum is continuous) their method 
provided the real part of the energy of the resonance only for a discrete set of parameters 
and this set could not be chosen a priori. Notwithstanding this, Ferron et ai, conjectured 
that there is a well defined function S{Er), the von Neumann entropy of the resonance state, 
which has a well defined value for every value of the real part of the energy of the resonance. 

In this work we will show, if A is the external parameter that drives the quantum dot 
through the ionization threshold, that the entropy S{Er{X)) is a smooth function of A. Also 
it is showed that the resonance state entropy calculated by Ferron et al. is correct near the 
ionization threshold. 

We have studied other quantities, besides the entropy, that are good witnesses of reso- 



nance presence. One of them is the Fidelity, which has been widely used (20|, I20l422l| in the 
detection of non-analytical behavior in the spectrum of quantum systems. The analysis of 
the Fidelity provides a method to obtain the real part of the resonance energy from vari- 
ational eigenstates. We introduce the Double Orthogonality function (DO) that measures 
changes in quantum states and detects the resonance region. The DO compares the ex- 
tended continuum states and the state near the resonance, also providing the real part of 
the resonance energy. 

The paper is organized as follows. In Section [III we present the model and briefly explain 
the technical details to obtain approximate eigenvalues, eigenfunctions, and the density of 
states for the problem. In Section IIIII the fidelity is used to obtain the real part of the 
resonance energy and the Double Orthogonality is introduced as an alternative method. In 
Section IIVI the linear entropy and the expectation value of the Coulombian repulsion are 
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studied using complex scaling methods. Finally, in Section |V] we discuss our results and 
present our conclusions. 



II. THE MODEL AND BASIC RESULTS 

There are many models of quantum dots, with different symmetries and interactions. In 
this work we consider a model with spherical symmetry, with two electrons interacting via 
the Coulomb repulsion. The main results should not be affected by the particular potential 
choice as it is already known that the near threshold behavior and other critical quantities 
(such as the critical exponents of the e nerg y and other observables) are mostly determined 



by the range of the involved potentials [2^. Therefore to model the dot potential we use a 
short-range potential suitable to apply the complex scaling method. After this considerations 
we propose the following Hamiltonian H for the system 

H = -^V^, - + V{n) + V{r,) + (1) 

2m ^ 2m ^ |r2 — ri| 

where V{r) = — (Vo/rQ) exp (— r/rg), the position operator of electron i = 1,2; Vq and Vq 
determine the range and depth of the dot potential. After re-scaling with tq, in atomic units 
the Hamiltonian of Eq. ([1]) can be written as 

H = -^Vl^ - - T/oe-- - K,e-- + ^^^-A-, (2) 

where A = rg. 

We choose the exponential binding potential to take advantage of its analytical properties. 
In particular this potential is well behaved and the energy of the resonance states can 
be calculated using complex scaling methods. So, besides its simplicity, the exponential 
potential allows us to obtain independently the energy of the resonance state and a check 
for our results. The threshold energy, e, of Hamiltonian Eq. (|2]), that is, the one body ground 



state energy can be calculated exactly 



25[ and is given by 



J2V2-S V^o) = 0, (3) 



where J^(x) is the Bessel function. 

The discrete spectrum and the resonance states of the model given by Eq. ([2]) can be 
obtained approximately using variational functions 17|, 23|. So, if \ipj{l,2)) are the 



exact eigenfunctions of the Hamiltonian, we look for variational approximations 



M 



1^,(1,2)) ^ |vl;,.(l,2)) = 5^cp)|<l>,) , c?) = (c(^-)). ; j = l,---,M. 



(4) 



i=l 



where the must be chosen adequately and M is the basis set size. 

Since we are interested in the behavior of the system near the ground-state ionization 
threshold, we choose as basis set s-wave singlets given by 



= |ni,n2;/) = (0„i(ri) 0„,(r2)), 3^o,o(^i. ^2) Xs , (5) 
where ^2 < ni, / < ?t,2, Xs is the singlet spinor, and the y^ QiVti, VL2) are given by 



rn=—l 



(6) 



i.e. they are eigenf unctions of the total angular momentum with zero eigenvalue and the 
are the spherical harmonics. Note also that q is a real function since it is symmetric 
in the particle index. The radial term {4>ni{fi)4'n2{f'2))s has the appropriate symmetry for a 
singlet state, 



where 



(0ni(ri)0n2(r2))^ 



{ni\n2) 



0ni('^l)0n2('^2) + 0ni (^2)0n2 (n) 



[2(l + (ni|n2)2)] 



2^ll/2 



(7) 



(8) 



and the (/)'s are chosen to satisfy {ni\n-\) = 1. The numerical results are obtained by taking 
the Slater type forms for the orbitals 



a 



2n+3 



1/2 



(9) 



^" ^ ' [(2n + 2)!_ 

where a is a non-linear parameter of the basis. It is clear that in terms of the functions 
defined in Eq. the variational eigenfunctions reads as 



^]/f)(l,2))=5^c«£)Kn2;/;a), 



(10) 



nin2l 

where ni > ^2 > / > 0, then the basis set size is given by 

A'' ni n2 

1 = , 

6 



N n\ n2 , 

^ = E E E 1 = + + 2)(^ + 3) 



(11) 



ni=0n2=0 1=0 



so we refer to the basis set size using both and M. In Eq. fllUI) we added a as a basis 
index to indicate that in general the variational eigenfunction is a-dependent. The matrix 
elements of the kinetic energy, the Coulombic repulsion between the electrons and other 
mathematical details involving the functions |ni,n2;/;«) are given in references 



27|. 



We only show here for completeness the matrix elements of the exponential potential in the 
basis of Eq. ©, 




FIG. 1: (color on-line) (a) the figure shows the behavior of the variational eigenvalues E^°^^ {\) 
(black lines) for = 14 and non- linear parameter a = 2. The red dashed line corresponds to the 
threshold energy e ~ —1.091. Note that the avoided crossings between the variational eigenvalues 
are fairly visible, (b) The figure shows the same variational eigenvalues that (a) (black lines) and 
the energy calculated using the complex scaling method (green line) for a parameter 6 = vr/lO. 

Resonance states have isolated complex eigenvalues, E^es = — F > 0, whose 

eigenfunctions are not square-integrable. These states are considered as quasi-bound states 
of energy E^ and inverse life time F. For the Hamiltonian Eq. ([2]), the resonance energies 
belong to the interval (£,0) [isj l. 

The resonance states can be analyzed using the spectrum obtained with a basis of 



functions (see 



19| and References therein). The levels above the threshold have several 



avoided crossings that "surround" the real part of the energy of the resonance state. The 
presence of a resonance can be made evident looking at the eigenvalues obtained numerically. 
Figure [1] shows a typical spectrum obtained from the variational method. This figure shows 
the behavior of the variational eigenvalues E'j"'' as functions of the parameter A. The results 
shown were obtained using = 14 and a = 2.0. The value of a was chosen in order to 
obtain the best approximation for the energy of the ground state in the region of A where 
it exists. The Figure shows clearly that for A < Xth — 1.54 there is only one bounded state. 
Above the threshold the variational approximation provides a finite number of solutions with 
energy below zero. Above the threshold there is not a clear cut criteria to choose the value 
of the variational parameter. However, it is possible to calculate -E^GV) calculating Ej"^ for 
many different values of the variational parameter (see Kar and Ho 28l|). 

Figure [2] (a) and (b) shows the numerical results for the first and second eigenvalues 
respectively, for different values of the variational parameter a. The figure also shows the 
behavior of the ground state (below the threshold) and the real part of the energy of the 
resonance calculated using complex scaling (above the threshold), this curve is used as a 
reference. The behavior of the smaller variational eigenvalue e["\X) is rather clear, below 
the threshold e["\X) is rather insensitive to the actual value of a, the differences between 
^i"='^(A) and ^J"='^(A) are smaller than the width of the lines shown in the figure. Above 
the threshold the behavior changes, the curve for a given value of a has two well defined 
regions, in each region Ei is basically a straight line. The two straight lines in each region 
has a different slope and the change in the slope is located around £',.(A). 

In the case of -Eg""* (A) there are three regions, in each one of them the curve for a given 
value of a is basically a straight line and the slope is different in each region. A feature 
that appears rather clearly is that, for fixed A , the density of levels for energy unit is not 
uniform, despite that the curves Ej°'^\X) are drawn for forty equally spaced a^'s between 
a = 2.0 and a = 6.0. This fact has been observed previously [33| and the density of states 
can be written in terms of two contributions, a localized one and an extended one. The 
localized density of states is attributed to the presence of the resonance state, conversely 
the extended density of states is attributed to the continuum of states between [e, 0). 



The localized density of states p{E) can be expressed as [28|, l33| 



PiE) 



dE{a) 
da 



-1 



(13) 
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FIG. 2: (color on-line) (a) The first variational state energy vs A, for different values of the 
variational parameter a. From bottom to top a increases its value from a = 2 (dashed blue line) 
to a = 6 (dashed orange line). The real part of the resonance eigenvalue obtained using complex 
scaling { 9 = vr/lO) is also shown (green line), (b) Same as (a) but for the second state energy. 



Since we are dealing with a variational approximation, we calculate 



("i+i) 



(A) -K|"'-^)(A) 



(14) 



Figure [3] shows the typical behavior of Pj{E) = p{Ej (X)) for several eigenenergies and 
A = 2.25. The real and imaginary parts of the resonance's energy, Er{X) and F respectively, 
can be obtained from p{E), see for example 28| and references there in. This method 
provides and independent way to obtain E^es, besides the method of complex scaling. 

The values of -E'r(A) and F(A) are obtained performing a nonlinear fitting of p{E), with 
a Lorentzian function, 



P{E) = Po 



A 



F/2 



(15) 



7r[(E-E,)2 + (F/2)2]- 
One of the drawbacks of this method results evident: for each A there are several Pj{E) 
(in fact one for each variational level), and since each Pi(E ) pr ovides a value for E^{X) and 
F-'(A) one has to choose which one is the best. Kar and Ho 28|] solve this problem fitting all 



the Pj{E) and keeping as the best values for Er{X) and F(A) the fitting parameters with the 




-0.8 -0.75 -0.7 



FIG. 3: (color on-line) The density of states p{E) for A = 2.25 and basis set size N = 14. The 
results were obtained using Eq. (jl3p and correspond to, from top to bottom, the second (black 
line), third (dashed red line), fourth (green line) and fifth (dashed blue line) levels. 



smaller value. At least for their data the best fitting (the smaller x^) usually corresponds 
to the larger n. This fact has a clear interpretation, if the numerical method approximates 
Er{X) with Ei"\X) a large n means that the numerical method is able to provide a large 
number of approximate levels, and so the continuum of states between (e, 0) is "better" 
approximated. 

It is worth to remark that the results obtained from the complex scaling method and 
from the density of states are in excellent agreement, see Table [H 



III. FIDELITY AND DOUBLE ORTHOGONALITY FUNCTIONS 



Since the work of Zanardi et al. 



20 



21| there has been a 



growing interest in the fidelity 



approach as a mean to study quantum phase transitions 20 , the information-theoretic 



differential geometry on quantum phase transitions (QPT's) 



XY model 



2l| or the disordered quantum 



22|. In all these cases the fidelity is used to detect the change of behavior of 



the states of a quantum system. For example, if A is the external parameter that drives a 
system through a QPT, the fidelity is the superposition J-" = (^(A — 6X), \E'(A + SX)), where 
\1/ is the ground state of the system. It has been shown that J-" is a good detector of critical 
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behavior in ordered [20|] and disordered systems |22[]. 

In the following we will show that the energy levels calculated using the variational 
approximation show critical behavior near the energy of the resonance, moreover the curve 
Er{X) can be obtained from the fidelity. 

Figure m shows the behavior of the function Qn = ^ — Fn-, where 
Fn = I (\1/„(A), \1/„,(A + 5A))p , and \l/„, is the n*'* eigenstate obtained with the variational 
approach. 

The behavior of Q is quite simple. The value of Q is very small, except near the avoided 
crossings where the value of Q increases rather steeply (at least for small n). This is so 
because near the avoided crossing the superposition |(\1/„(A), \I^ri(A + 5A))p — )■ 0. Actually, 
|(\1/„(A), \1/„(A + 5A))p — )■ near points such that -E'n"^(A) has non analytical behavior. Is 



for this reason that the fidelity is a good detector of quantum phase transitions [20|-|22|. In 
a first order QPT the energy of the ground state is non analytical, and in a second order 
QPT the gap in the avoided crossing between the ground state and the first excited state 
goes to zero in the thermodynamic limit. 

The previous argument supports why Qi has only one peak, while all the others functions 
Qn have two peaks, the number of peaks is the number of avoided crossings of each level. 
However, since the resonance state lies somewhere between the avoided crossings it is natural 
to ask what feature of the fidelity signals the presence of the resonance. For a given level n 
the value of the energy is fixed, so we must pick a distinctive feature of Qn that is present 
for some A{, such that Er{\{) — -E'„(A{) (from here we will use En or Elf* interchangeably). 
It results to be, that A^ is the value of A such that Qn attains its minimum between its two 
peaks. Figure m shows the points En{X-l). In Table [I] we tabulate the real part of the energy 
calculated using DO, complex scaling, fidelity and Density of States, for the five values of A{ 
shown in Figure HI The numerical values obtained using the fidelity and Density of States 
methods are identical up to five figures. The Relative Error between the energies obtained 
is less than 0.25% 

The idea of detecting the resonance state energy with functions depending on the inner 
product could be taken a step further. To this end we consider the functions 

DOn{\) = |(vl>„(Ai), vl>,(A))|2 + |(vl^„(A^),vI/„(A))|2,for Ai < A < Ar (16) 

where A^ and Xr are two given coupling values. It is clear from the definition that 
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-0.75 - 



-1.25 




FIG. 4: (color on-line) The upper panel shows the behavior of for n = 1, . . . ,7. Each function 
Gn has two peaks, except for n = 1. Since one of the peaks of Gn coincides with one of the peaks 
of Qn+i, only one peak for each level n is apparent. From left to right the visible line at each peak 
correspond to n = 2, . . . , 7 (cyan, red, yellow, blue, grey, brown lines, respectively). The n = 1 line 
has no visible peak (black) . The lower panel shows the variational eigenlevels n = 1 , . . . , 7 (with 
the same color convention used in the upper panel), and Er{X) (green dashed line). The black 
vertical dashed lines connecting both panels show the value of A where each Gn has its minimum. 
An- The red dots in the lower panel correspond to En{Xn). 



< -DO„(A) < 2. If there are not resonances between Xl and Xr, the wave function is 
roughly independent of A so DOn{X) ~ 2. However, the scenario is different when a resonance 
is present between Xl and Xr. In this case, the avoided crossings for a given state are lo- 
cated approximately at A^^ and A^^, where L{R) stands for the leftmost (rightmost) avoided 
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TABLE I: Resonance Energy obtained by four different methods. Basis size is N=14 









DO 


Complex 
Scaling 


Fidelity and Density of States 






n = 2 


n = 2> 


n = 4 n = 5 


n = 6 


1.755 (n= 


=2) 


E 


-0.99075 


-0.99098 


-0.99011 


— 


— — 


— 






Q! 


2.0 


2.0 


1.560 


— 


— — 


— 


1.8625 (i) 


=3) 


E 


-0.93427 


-0.93452 


-0.93434 


-0.93383 


-0.93303 — 


— 






OL 


2.0 


2.0 


2.448 


1.787 


1.414 — 


— 


2.02 {n= 


=4) 


E 


-0.85498 


-0.85556 


-0.85581 


-0.85564 


-0.85531 -0.85486 








a 


2.0 


2.0 


3.339 


2.435 


1.906 1.519 




2.255 (n= 


=5) 


E 


-0.74329 


-0.74538 


-0.74518 


-0.74527 


-0.74519 -0.74521 


-0.74514 






a 


2.0 


2.0 


4.262 


3.098 


2.414 1.936 


1.574 


2.61 (n= 


=6) 


E 


-0.58276 


-0.59077 


-0.58825 


-0.58910 


-0.58942 -0.58965 


-0.58979 






a 


2.0 


2.0 


5.248 


3.799 


3.799 2.373 


1.936 



crossing. Requesting that < A^'' < A^" < Ar, it follows that (^n(AL)|^„(Aij)) ~ 0. With 
this prescription, the DOn functions are rather independent of the actual values chosen for 
\l and \r. For a given n, DOn{X) measures how much the state \l/n(A) differs from the 
extended states \1/„(Ar) and '^ni.^Lj- 

We look for the states with minimum DOn-, in the same fashion as we did with the 
fidelity, we can obtain values Er{X^o) — ^ni^Do)' where A^^o is defined by DOn{XDo) = 
minA-DOn(A). Figure [5] shows the behavior of DOn obtained for the same parameters that 
the ones used in FigureHl and we compare the values of £J„(A^q) with energy values obtained 
using complex scaling methods in TablelB The curves in FigureOshow that outside (A^,, Ar) 
the states change very little when A changes and DOn — 1- Inside the resonance region, 
{X'l', A^*"), the functions DOn change abruptly. The width in A in which a given DOn changes 
abruptly apparently depends on the width of the resonance, F, but so far we have not been 
able to relate both quantities. 

From Table [T] and Figure \5\ it is rather clear that despite that the fidelity and the DOn 
provide approximate values for Er{X) for different sets of A's, both sets belongs to the "same" 
curve, i.e. the same curve considering the numerical inaccuracies. Both methods would give 
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FIG. 5: (color on-line) The lower panel shows the variational energy levels £'„(A), from bottom 
to top, for n = 1,...,7 (the black, green, red, yellow, blue, grey and brown continuous lines, 
respectively); Er{X) (the dashed dark green line); En{y}jo) (blue squares) and En{\i) (red dots). 
The upper panel shows the behavior of DOn vs A for n = 2, ... ,7. The color convention for the 
DOn is the same used in the lower panel. The black dot-dashed vertical lines show the location of 
the points X^q- 

the same results when |A^*' — A^''] — 0, but for finite the fidelity measures how fast the 
state changes when A — A -|- AA and the DO measures how much a state differs from the 
extended states located at both sides of the resonance state. 
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IV. THE ENTROPY 

If jf^'^ is the reduced density operator for one electron [19| , then the von Neumann entropy 
S is given by 

s = -ti{r''^og,r^), (17) 

and the hnear entropy is given by 11 1 



5u,= l-tr[(p-'i)2], (18) 
where the reduced density operator is 

r\v,,v\) = ti^\^){m\ , (19) 

here the trace is taken over one electron, and |\E') is the total two-electron wave function. 
Both entropies, Eqs. ( 1T7|) and ( IT8|) . can be used to analyze how much entanglement has a 
given state. One can choose between one entropy or the other out of convenience. In this 
paper we will use the linear entropy. For a discussion about the similarities between the two 



entropies see Reference [ll[ and references therein. 

As the two electron wave function is not available we instead use the variational approx- 
imation Eq. (fTOj) . As has been noted in previous works (see jsj and References therein), 
when the total wave function factorizes in spatial and spinorial components it is possible to 
single out both contributions, then the analysis of the behavior of the entropy is reduced to 
the analysis of the behavior of the spatial part S*, since the spinorial contribution is constant. 
In this case, if v?(ri, r2) is the two electron wave function and p'''^'^(ri, r'^) is given by 

r\v,, v\) = I ^*(ri, T2Mr[, rs) dr^, (20) 
then the linear entropy Su^ can be calculated as 

5h„ = 1-$^A2, (21) 

i 

where the A, are the eigenvalues of p'^'^'^ and are given by 

I p'^'iru r[)<P,{r[) dr[ = A,0.(ri) . (22) 

Figure [6] shows the behavior of the linear entropy for several variational levels. The 
meaning of each curve has been extensively discussed in Reference 19|]. We include a brief 
discussion here for completeness. 
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FIG. 6: The figure shows the behavior of S'iin(^'^"^^ 



.J where the are the variational eigenstates 
corresponding to the first seven energy levels shown in Figure [1] for = 14 and a = 2.0. All the 
curves 5'iin(^j"^) , except for the corresponding to 5iin(^i"^), have a single minimum located at 



Aj"', i.e. ^(Af)) = minASiin(*f^(A)) . If z < j then Af < Af. 



When the two-electron quantum dot loses an electron the state of the system can be 
described as one electron bounded to the dot potential, and one unbounded electron at 
infinity, as a consequence the spatial wave function can be written as a symmetrized product 
of one electron wave functions so ^lin = S'c = 1/2. Therefore if only bound and continuum 
states are considered the entropy has a discontinuity when A crosses the threshold value 
\th- The picture changes significantly when resonance states are considered. The resonance 
state keeps its two electrons "bounded" before the ionization for a finite time given by the 
inverse of the imaginary part of the energy. Of course the life time of a bounded state is 



infinite while the life time of a resonance state is finite. In reference [19| is suggested that it is 
possible to construct a smooth function S{Er{\)) that "interpolates" between the minima of 
the functions j) shown in Figure El This assumption was justified by similar arguments 
that the used in the present work, i.e. if we call Af to the value of A where 5'(\l/„,] gets its 



19| used 



minimum then i?„(Af) follows approximately the curve Er{\). As Ferron et al. 
only one variational parameter a, it seemed natural to pick the minimum value of S{^n) as 
the feature that signaled the presence of the resonance state. 
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Until now we have exploited the fact that Er{X), at a given A can be approximated 
by variational eigenvalues corresponding to different values of the variational parameter, 
say Er{X) — Ei^\\) ^ E^^\X) (the superscript a is made evident to remark that the 
eigenvalues correspond to different variational parameters a and a' ). There is no problem 
in approximating an exact eigenvalue with different variational eigenvalues. But, from the 
point of view of the entropy, there is a problem since, in general, S{^^\\)) is not close 
to S'(\E'n" ^(A)). Moreover, as has been stressed in Reference [2I, a given numerical method 
could be useful to accurately calculate the spectrum of a quantum system, but hopelessly 
inaccurate to calculate the entanglement. In few body systems there is evidence that there 
is a strong correlation between the entanglement and the Coulombian repulsion between 

n n n n 

the components of the system [5|, l9|, lUl, l26| . Because of this correlation we will carefully 
investigate the behavior of the Coulombian repulsion between the electrons in our model. 

For the Hamiltonian Eq. ([2]), and ^ G £^ an eigenvector of H with eigenvalue E, the 
Hellman-Feynman theorem gives that 

n 771 1 

(23) 

We use both sides of Eq. (1251) to analyze how the variational approximation works with 
expectation values of observables that are not the Hamiltonian. The r.h.s of Eq. fl23l) is 
well defined if we use functions as the approximate variational eigenfunctions. 

To evaluate the l.h.s of Eq. fl23|) we take advantage that we have found, independently, the 
real part of the resonance eigenvalue, Er{X), using complex scaling methods. Figure [7] shows 
the behavior of and the Coulombian repulsion between the two electrons, (^)n; where 
(...)„ stands for the expectation value calculated with \E'i"'*. The behavior of (^)n is quite 
simple to analyze, where the linear entropy of \E'i°'' has a valley the expectation value (^)n 
has a peak. Where the expectation value {:;^)n has its maximum the corresponding linear 
entropy has its minimum. The inverse behavior showed by the entropy and the Coulombian 
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26| 



repulsion has been observed previously 

For a given variational parameter a, and for small n, (^)n has its maximum very close to 
the curve Besides, the shape of both curves near the maximum of (^)n is very similar, 
in this sense our variational approach gives a good approximation not only for Er{X) but 
for its derivative too. 

For larger values of n the maximum of (;:^)n gets apart from the curve of and the 
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FIG. 7: The figure shows the expectation values of the Coulombian repulsion for the variational 



states ^ 



(a) 



n 



1, ... ,8 with iV = 14 and a = 2.0. Also is showed the curve obtained from 



the complex Energy of the complex scaling method. 



shape of the curves near this maximum is quite different. We proceed as before and changing 
a we obtain a good approximation for up to a certain value A^ep- For any A smaller 
than Xrep, there is a pair n, a such that (^)„,a is locally close to ^ and the slope of both 
curves is (up to numerical errors) the same, see figure [HI 

Apparently there is no way to push further the variational method, at least keeping the 
same basis set, in order to obtain a better approximation than the depicted in Figure [HI The 
difficulty seems to be more deep than just a limitation of the variational method used until 
this point. We can clarify this subject using the properties of the complex scaling method. 
Let us call 0^ the eigenvector such that 



(24) 



where H{6) is the Hamiltonian obtained from the complex scaling transformation 29|, and 



6 is the angle of "rotation". The eigenvector 0^ depends on 6, but 



eigenvalue Eres does not depend on 6. As pointed by Moiseyev 



br 6 large enough the 



29| . the real part of the 



expectation value of a complex scaled observable is the physical measurable quantity, while 
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FIG. 8: (color-online) The figure shows the expectation value {l/ri2)^'^ vs. A for a basis size 
= 14 and a = 2, ... , 3.5 in 0.1 steps and for a = 4, ... , 5.5 in 0.5 steps (solid black lines). The 
real (cyan dotted) and imaginary (orange line) parts of {'i/ri2)e {0 = tt/IO), and the derivative of 
the real part of the complex-scaled energy are also shown. 



the imaginary part gives the uncertainty of measuring the real part. Moreover, the physical 
measurable quantity must be 9 independent as is, for example, the eigenvalue E^es- 
The eigenvector 0^ can be normalized using that 



(25) 



Since 0^ is normalized, we get that 

BE 



-ie 



1 

ri2 



dX ^ ?^i2 

in this last equation we have used that, under the complex scaling transformation. 



(26) 



1 e 
> - 



Ad 



and defined the quantity ( ^ 



for Gamow states 



30|. 



(27) 

ri2 ri2 

This generalized Hellman-Feynman theorem is also valid 



Figure [8] shows the behavior of the expectation value in the ( — / as a function of A. 



It is clear that the real part of the expectation value coincides with More 
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interestingly, A^ep is where the imaginary part of 




became noticeable. From this fact 



we conclude that it is not possible to adequately approximate the Coulombian repulsion of 
a resonance state, or its entropy, only with real variational functions despite its success 
when dealing with the resonance state spectrum. 

We define the complex scaled density operator of the resonance state by 

/ = |/>((0Y|, (28) 

and the complex linear entropy 

S' = l-ir{plJ\ (29) 

where 

Pld = tr2p', (30) 

and (f)^ is the eigenvector of Eq. f l2^ . This definition is motivated by the fact that the den- 
sity operator should be the projector onto the space spanned by As the normalization 
Eq. (l25i) requires the bra to be conjugated, then is the adequate projector to use. 

Because of the normalization, Eq. f l25|) . we have that trp^ = trpf^^ = 1 despite that both 
density operators have complex eigenvalues. 

Figure M shows that up to certain value of A the real part of follows closely a enve- 
lope containing the minima of the functions S'(\E'„). However, for A large enough, gets 
apart from the functions S'(\E'n)- It is worth to mention that for 9 large enough does 
not depend on 6. On the other hand, far away from the threshold, the complex scaling 
requires larger values of 6 to isolate the resonance state eigenenergy, but in this regime the 
method becomes unstable. Because of the numerical evidence, near the threshold the en- 
tropy calculated by Ferron et al. is basically correct, but for larger values of A the amount 
of entanglement of the resonance state should be characterized by 5*^ and not by any of the 



V. SUMMARY AND CONCLUSIONS 



We have presented numerical calculations about the behavior of the fidelity and the double 
orthogonality functions Z^O„(A). The numerical results show that it is possible to obtain 
Er{X) with great accuracy, for selected values of A, without employing any stabilization 



19 




FIG. 9: (color-online) Figure (a) shows the linear entropy for the same values as in figure [6] (magenta 
dashed lines). Also showed is the real part of the complex linear entropy for several values of the 
complex rotation angle ^ = f > fgi ^) g^; (black empty diamonds, red dots, green squares, blue 
triangle, yellow empty dots). Figure (b) shows the imaginary part of the complex linear entropy 
for the same values as in (a). 

method. These two methods to find Er{X) do not depend on particular assumptions about 
the model or the variational method used to find approximate eigenfunctions above the 
threshold, their success depends on the ability of the approximate eigenstates to detect the 
non-analytical changes in the spectrum. 

The fidelity has been extensively used to detect quantum phase transitions in spin sys- 
tems [2^, the behavior of quasi-integrable systems 31], thermal phase transitions 32], etc. 
To the best of our knowledge this work is the first attempt to apply the concept of fi- 
delity to resonance states and to the characterization of spectral properties of a system 
with non-normalizable eigenstates. Besides, it is remarkable that the fidelity and the dou- 
ble orthogonality give the real part of the resonance eigenvalue using only real variational 
functions. This energy as a function of A is obtained by moving the nonlinear parameter a 
but without fitting as needed by standard stabilization methods. Moreover, as shown by the 
tabulated values in Table [T] the fidelity provides Er{X) as accurately as the density of states 
method, with considerable less numerical effort. 

We proposed a definition of the resonance entropy based on a complex scaled extension of 
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the usual definition. The extension imphes that the reduced density operator is not hermitian 
and has complex eigenvalues, resulting in a complex entropy. The real and imaginary parts 
of the complex entropy are 9 independent, as should be expected for the expectation value 



of an observable [29|. This independence gives support to the interpretation of the real part 
of the entropy as the amount of entanglement of the resonance state. 

Other kinds of resonances, as those that arise from perturbation of bound states embedded 
in the continuum, could be studied applying the same quantum information methods used 
in this paper. Work is in progress in this direction. 
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